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Abstract 

Pairing heaps are shown to have constant amortized time Insert and Meld, thus 
showing that pairing heaps have the same amortized runtimes as Fibonacci heaps for 
all operations but Decrease-key. 



1 Introduction 
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Pairing heaps were introduced in [10] as a priority queue data structure modeled after splay 
As in splay trees, pairing heaps do not augment the nodes of the heap with any 



trees 

information, and use simple local restructuring heuristics to perform all operations. They are 
easy to code and empirically perform well [l9] . 

While the study of splay trees has been focused on their ability quickly execute various 
distributions of operations [l}]6| [l3|[l5|[l8|[22] , the study of pairing heaps remains stuck at a 
much earlier stage — tight amortized bounds on the runtimes of all operations in terms of 
the heap size remain unknown. There is also small body of work to study how pairing heaps 
work on particular distributions of operations [7,12,16 



On the practical side, pairing heaps are in use. They are covered in some elementary and 

Additionally, one can find on the web code 



intermediate level data structures texts 23 -26 



that has been developed to implement pairing heaps in various languages. Pairing heaps were 
part of the pre-STL GNU C++ Library, and thus were distributed widely. 

The theoretically leading non-self adjusting priority queue is the Fibonacci heap 11 , and 
pairing heaps are designed to support the same set of operations at Fibonacci heaps: 

• /i=Make-Heap(): Returns an identifier h of a new empty heap. 

• x=Extract-MIN(/i): Removes and returns x, the minimum element in heap h. 
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*Part of this work appeared in preliminary form in |12 and 
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Figure 1: Summary of previous results, in chronological order of first appearance. All 
previous results support 0(1) amortized Make- heap and O(logn) amortized Extract-min. 
All results are amortized, so the runtimes of the individual operations can not be mixed 
among the different results. The result of Fredman's [9] gives an amortized lower bound for 
Decrease-key, given amortized upper bounds for other operations. The result of Stasko 
and Vitter is for a variant of pairing heaps and does not allow the Decrease-KEY operation. 



• p=Insert(/i, x): Inserts x into heap h, and returns an identifier p that can be used to 
manipulate x in the future. 

• Delete(/;,,p): Removes the item in h identified by p. 

• Decrease-key(/i,p, A): Decreases the key value of the item in h identified by p by a 
nonnegative amount A. 

• /;,=Meld(/ii, /;,2): Combines the contents of heaps hi and h2 into a new heap with 
returned identifier h. The identifiers to hi and /12 are no longer valid. 

• x=FIND-MIN(/;,): Returns the minimum element of the heap h. 

In the original pairing heap paper |10|, all operations except the constant-time Make- 
heap and FiND-MiN were shown to take 0(log?T,) time. It was conjectured in |lO] and 



empirical evidence was presented by Stasko and Vitter 19 that pairing heaps share the 
same amortized cost per operation as Fibonacci heaps, which have 0(1) decrease-key, 
insert and MELD operations. However, this possibility was eliminated when it was shown 
by Fredman [9j that the amortized cost of decrease-key can not be below O (log logn). 
Recently, Pettie has produced a new analysis that focuses on the decrease-key operation, 
where he proves a 0(2^^^°^^°^") amortized bound on insert, meld, and decrease-key 
while retaining a O(logn) upper bound on extract-min. 

In this work we present a new analysis of pairing heaps that proves, with the exception of 
DECREASE-KEY operations, pairing heaps share the same asymptotic runtime per operation 
as Fibonacci heaps. Specifically, we show the amortized cost of extract-min, delete, and 
DECREASE-KEY is O(logn), the amortized cost of make-heap and insert is 0(1), and the 
amortized cost of meld is 0. Thus, compared to the original analysis in |10| the amortized 
upper bound of O(logn) for the insert and meld operations, is improved to 0(1) for insert 
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and cj^for MELD. Compared to the analysis of Pettie, the 0(2^^'°s^°g") bounds for insert 
and MELD are improved to a constant, while Pettie's 0(2^^^°^'°^"-) bound for decrease-key 
is tighter than our 0(logr;,) bound. 

It should be noted that Stasko and Vitter in introduced a variant of pairing heaps, 
the auxiliary twopass method, and proved that this structure supported constant amortized 
time Insert. However, their analysis explicitly forbade the Decrease-key operation. 
Elmasry invented a different variant of pairing heaps that obtains amortized 0(1) Insert 
and amortized O(loglogn) Decrease-KEY; however Fredman's Q(loglogn) amortized lower 
bound on Decrease-key does not apply to Elmasry's variant as they do not confirm to his 
model of generalized pairing heaps. 

Here is is shown that pairing heaps are shown to have constant amortized time Insert 
and Meld, thus showing that pairing heaps have the same amortized runtimes as Fibonacci 
heaps for all operations but Decrease-KEY. The Decrease-KEY is allowed at an amortized 
cost of 0(logr;,). 

2 Pairing Heaps 

A pairing heap is a heap ordered general tree. Here a min-heap is assumed. The basic 
operation on a pairing heap is the pairing operation, which combines two pairing heaps into 
one by attaching the root with the larger key value to the other root as its leftmost child. 
Priority queue operations are implemented in a pairing heap as follows: Make-heap creates 
a new single node heap. Find-min returns the data in the root of the heap. Meld pairs 
the roots of the two heaps. Insert pairs the new node with the root of heap it is being 
inserted into. Decrease-key breaks off the node and its induced subtree from the heap (if 
the node is not the root), decreases the key value, and then pairs it with the root of the heap. 
Delete breaks off the node to be deleted and its subtree, performs an an extract-min on 
the subtree, and pairs the resultant tree to the root of the heap. Extract-min is the only 
non-trivial operation. An Extract-min removes and returns the root, and then, in pairs, 
pairs the remaining trees in the resultant forest. Then, the remaining trees from right to 
left are incrementally paired. All pairing heap operations take constant actual time, except 
EXTRACT-MIN and DELETE, which take time linear in the number of children of the node to 
be removed. For the purposes of implementation, pairing heaps are stored as a binary tree 
using the leftmost child, right sibling correspondence. Unless otherwise stated, the standard 
tree terminology will refer to the general tree representation. 

^Meld in Fibonacci heaps is typically stated as taking 0(1) amortized time. However, since meld 
operations must be dominated by make-heap operations, valid meld operations can never asymptotically 
change the runtime of any sequence, and thus take amortized time. 
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3 Main Result 



3.1 Overview 

We claim that in a pairing heap the amortized runtime of find-min, make- heap, and insert 
is 0(1), MELD is and decrease-key, delete and extract-min is 0(logr;,). We adopt 
the convention that log refers to the binary logarithm. 

Let X = xi,X2, . . . ,Xm be a sequence of operations to be executed an initially empty 
collection of heaps. The remainder of the notation used is implicitly parameterized by a fixed, 
but arbitrary, X. Let be the actual time to execute operation Xj. Let be the size of the 
heap that Xi acted upon after the execution of Xj. We define the actual cost of the operation 
to be the number of pairings performed plus 1. This cost measure is chosen as the number of 
parings clearly dominates the asymptotic runtime of each non-constant cost operation. Let Oj 
be the amortized cost of operation Xj. Specifically, let cij = if Xj is a meld, a, = 0(1) if Xj is 
an INSERT, find-min, or a make-heap, and let = O(lognj) if x, is a delete, decrease-key 
or extract-min. 

Given this notation, the main theorem can be stated. The time to execute any sequence 
on an initially empty pairing heap is bounded but the sum of the previously stated amortized 
times for each operation. Formally, 

Theorem 1. YlT=i(^i ^ YlT=i^i 

The potential method is used to prove this theorem. We define below a potential function 
$ = $0; ^1; • • • ; that Is slmply a sequence of real numbers. The potential method may be 
summarized in the following lemma: 

Lemma 2. // there exists a sequence $ such that for all i, 1 < i < m, cij > + $j — 
and $m - $0 > then Yl'iLi > YllLi ^i- 

This lemma, which summarizes the potential method, may be proved by simple algebraic 



manipulation. More details on the potential method may be found in 21 



The proof of the theorem 1 proceeds as follows. Section 3.2 is devoted to defining the 
potential function $. This potential function is complex, and was derived using trial- and- 
error. We then prove, in a sequence of lemmas, that > + $j — for each type of 
operation Xj. Finally, we state in a lemma that — $o > 0. These lemmas, according to 
Lemma [2] are sufficient to prove Theorem [TJ Analysis of Insert is not presented separately, 
as Insert is just a Make-heap followed by a Meld. Similarly, the analysis of Delete is 
not presented separately, as Delete its analysis follows directly from that of Decrease-key 
and Extract-min. 



3.2 The potential function 

For the analysis, a color, black or white, is assigned to every node, and a weight is assigned to 
those nodes colored white. A node is black if it will remain in the forest of heaps at the end 
of execution of sequence X, and white otherwise. The color of a node can never change. We 
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say that a white node has a weight of heavy if the number of white nodes in its left subtree in 

the binary representation is greater than or equal to the number of white nodes in its right 
subtree. Roots are always heavy by this definition and every node in a heap of size n can 
have a maximum of log n heavy children. We say that the weight of a white node that is not 
heavy is light. 

We say a node has been captured if its parent is black. A captured node must have a 
DECREASE-KEY performed on it later in the execution sequence before it is involved in any 
pairings. White captured nodes must have DECREASE-KEY or a DELETE performed on them 
later in the execution sequence. 

The potential of a white node is the sum of four components: rank potential, weight 
potential, triple white potential, and capture potential. Let s{x) be the number of white 
nodes in the induced subtree of x in the binary representation. The rank potential of a white 
node X, r{x) is 91ogs(a;). If node is white and has right and left siblings that are also white, 
then the node is referred to as a triple white and has a triple white potential of 0, else it has 
6 units of triple white potential. White heavy nodes have a weight potential of 0, and white 
light nodes have 6 units of weight potential. We assign captured nodes a capture potential of 
and non-captured nodes 6 units of capture potential. Black nodes are defined to have 
rank potential, weight potential and triple-white potential. Thus, the potential of a black 
node is only comprised of its capture potential. We also assign each heap a heap potential 
which is 8 — 18 Ya=i ^ogi where n is the number of white nodes in the heap. The potential of 
a forest of heaps, $j, is the sum of the potentials of the nodes in the heaps and the potentials 
of the heaps as a function of the state of the pairing heap after the execution of Xi. 

It must be noted that the potential function exists for the sole purpose of analysis. To 
implement a pairing heap neither the potential function, nor its constituents such as node 
color, need to be stored. 

The amortized cost of each operation is now calculated using this potential function. For 
each operation it is proved that > a, -|- $i — or in English, the amortized cost of an 
operation is at most the actual cost plus the change in potential. In order to analyze the 
the change in potential, we analyze the change in each of the components that are linearly 
combined to make the potential function. For each operation, only a limited number of clearly 
defined nodes may have their node potential change. These are the nodes that have had their 
parents change, or white nodes that have had a neighboring sibling or white descendant in 
the binary representation change. 

3.3 Make-heap 

Lemma 3. If Xi is a make-heap, ai + ^i- < 21. 
Proof. 

Actual cost: 1. No pairings are performed. 

In a MAKE- HEAP, the only change in potential is caused by the introduction of the new 
one-node heap. The potential of all existing nodes and heaps is unchanged. 
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Change in rank potential: 0. The newly created root, if white, has 1 node in its induced 
subtree, and thus has 9 log 1 = units of rank potential. If the newly created root is 
black, black nodes by definition have zero rank potential. 

Change in weight potential: 0. If the newly created node is white, it will be heavy and 
will have units of weight potential. If the newly created node is black, it does not 
have weight potential. 

Change in triple white potential: 6. The new node is not a triple-white. 

Change in capture potential: 6. The new node is not captured. 

Change in heap potential: 8 The new heap has a potential of 8. 

Summing the actual cost and the change in potential yields the amortized cost: Oj + — 
^i-i = 21 

□ 

3.4 Meld 

Lemma 4. If Xi is a meld, + $j — < 0. 

Proof. The only nodes in a meld that can change potential are the two old roots, and the 
leftmost child of the new root. We use a and b to denote the number of white nodes in the 
heaps being melded. 

Actual cost: 2, as 1 pairing is performed. 

Change in weight potential: < 6. Only new root can have its weight potential change. 

Change in triple white potential: < 0. Only the leftmost child of the root which will 
become the root of the resultant tree can change its triple white potential. However this 
change can only be negative, since this node is not a triple white before the beginning 
of this operation. 

Change in capture potential: < 0. No gain is possible since no node escapes capture in 

a MELD. 

The remainder of the analysis breaks into two cases. 

Case 1: At least one of the two heaps being melded only contains black nodes. 

Change in rank potential: 0. No white nodes have any change in their white 
descendants. 

Change in heap potential: -8. The heap with all black nodes (or one of them 
if both have only black nodes) has a heap potential of 8, while the resultant 
heap has the same potential as the other heap. 
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Case 2: Each of the heaps being melded contain some white nodes. 

Change in rank potential: < 181og(a + b). The only two nodes that could 
possibly change rank potential would be the two roots, and their rank potential 
could rise to be at most 9 log(a + b) each. 

Change in heap potential: < — 181og(a + b) — 8. 

The change in heap potential is given by the expression: 
8-18 E -i' log z - 8 + 18 Er=i log ^ - 8 + 18 E ti log t 
= -8 - 18 Er=f log i + 18 Eti log i + 18 E -=2 log i 

< -8 - 18 E ■=? log i + 18 Er=f log i 

< -8 - 181og(a + 6) 

Summing the actual cost and the change in potential for both Case 1 and Case 2 yields 
the same amortized cost: + — < 0. □ 

3.5 Decrease- key 

Lemma 5. If Xi is a decrease-key, + - < 26 + 15 log n^. 

Actual cost: 2, since one pairing is performed. 

Rank Potential: < Qlogn^. The node on which the decrease-key is performed could 
gain as most as 9 log rii in rank potential. 

Weight potential: < Glogrij + 12. The node the decrease key is performed on can gain 
6 in weight potential. Also, on the path from the node on which the decrease-key is to 
be performed to the root, the removal of the node and its subtree may cause some nodes 
to change their status from light to heavy or vice versa. Only changing from heavy to light 
causes a potential gain, and this gain of 6 can only happen in at most log + 1 nodes. 

Capture potential: < 6. There can only be a change of 6 in capture potential if the node 
which the decrease-key is performed on was captured. 

Triple white: < 6. Among the node on which the DECREASE-KEY is performed, and its 
two former siblings to the left and right, a total of 6 units of triple white potential can be 
gained. 

Summing the actual cost and the change in potential yields the amortized cost: Oj -|- $i — 
^i-i < 26 + 15 log n. 

3.6 Extract-min 

Lemma 6. Amortized cost of extract-min is 391og(ni + 1) + 89; 

Proof. Actual cost: If there are c children of the root, the actual cost is c, since c — 1 pairings 
are performed. 

As the analysis of extract-min is long, it is split into several parts. The first part 

examines the change of heap potential. Next, we note that the only nodes that can have their 
node potential change are the old root which is removed, the children of the old root, and the 
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grandchildren of the old root. The second part bounds the changes of node potential of the 
old root. The third part bounds the changes of triple-white potential of the grandchildren of 
the root, since this is the only type of node potential change possible in these nodes. The 
analysis of the node potential change in the children of the old root is the most complicated 
step, and it is presented in two parts: global and local. In the global part we analyze some 
potential changes over the whole structure, while in the local section we analyze how blocks 
of six children of the root are affected by the remaining potential changes. In the global 
section we analyze changes in rank potential, gains in weight potential, and gains of capture 
potential. In the local section we analyze losses of weight potential, and changes of triple 
white potential and losses of capture potential. 

Heap potential: 181og(nj + 1). This is caused by the heap's size being reduced from n + 1 
to n. 

Node potential of root: < 0. The removal of the root itself causes no potential gain, since 
it has nonnegative potential 

Triple white potential of grandchildren of the old root: < 0. Some grandchildren of 
the old root that are white, have white right siblings and no left siblings may, through 
a pairing that their parent is involved in, acquire a left white sibling and become a 
triple-white. This can only cause a loss of potential. 

Global Change of rank potential: < 27 logrii — lSw. Given that there are white- white 
pairings in the first pairing pass, the first pairing pass causes a rank potential 
gain of at most ISlognj — 18w; the derivation of this is a variation of the original 



pairing heap analysis of |10|. This is presented in Section 3.8 The second pairing 
pass causes a rank potential gain of at most 9 log n. 

Gains in weight potential: < Glogrij + G. There are at most logra + l heavy children 
of the root, and so the potential gain caused by heavy nodes becoming light is at 
most 6 log 77, j -|- 6. 

Gains of capture potential: 0. The extract-min operation can cause no increase 
in capture potential. An increase in the capture potential can only happen when 
a previously captured node becomes uncaptured. However, since the root is white, 
none of the children of the root can have black parents, and thus none were 
captured. 

Local: In order to analyze other changes in potential (changes in triple white potential, 
losses of weight potential caused by a node becoming heavy, and changes in black nodes' 
potential) we break the children of the root into blocks of six nodes, excluding the 
rightmost two nodes. At most 7 nodes can not be included in this analysis, and they 
could incur a potential gain of up to 12 each. In analyzing these specific potential 
changes in each block of six, there are six cases. 

Since for each case we are only considering losses in weight potential and capture 
potential, in some parts we do not explicitly state them and simply assume they are 
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nonpositive. Triple white potential must be considered in each part because we are 
considering gains, as well as losses, however we show gains are limited to Case 1. 

Case 1: There is at least one white- white pairing in the first pairing pass. 

Triple while: < 36. The only gains in potential are the possible gain of 6 units 
of triple white potential for each white involved in a white-white pairing. This is 
the only case where a gain in potential, among the components of the potential 
function under consideration is possible. A gain in triple white potential can only 
occur when a child of the former root is a triple white, and thus has white left 
and right siblings. However, this guarantees that any triple white be involved in 
a white-white pairing and this fall into this case. Thus, increases in triple white 
potential can not occur in the following cases and will not be considered. 

Case 2: There are no white- white pairings and at least one black-black pairing in the 
first pairing pass. 

Capture potential: < —6. The black-black pairing(s) causes a loss of at least 6 
units of potential. 

Case 3: All are black-white pairings in the first pairing pass, and at least one of the 
three white nodes is captured. 

Capture potential: < —6. The capturing of the node(s) causes a capture potential 
loss of at least 6. 

Ccise 4: All are black-white pairings in the first pairing pass, but all three nodes that 
participate in the second pairing pass lose. 

Triple white potential: < —6. Having all three loose the pairings in the second 
pairing pass causes a loss of potential of 6, due to the change of status of the 
middle white node to a triple white. 

Case 5: All are black- white pairings in the first pairing pass, and at least one of the 
three nodes that participate in the second pairing pass wins, and at least one of 
the white nodes is hght. 

Weight potential: < —6 The hght node becomes heavy, as all nodes previously on 
its right are now in its subtree. Additional nodes that were to the node's left may 
also be added to its subtree, but this just makes it more heavy. This causes a loss 
of 6 units of heavy potential. 

Case 6: All are black-white pairings in the first pairing pass, and at least one of the 
nodes that participates in the second pairing pass wins, and all of the whites that 
win in the second pairing pass are heavy. 

We may assume no potential gain. This case can only happen logn^ times, because 
there are at most logn^ heavy children of the root. 

Case 1 causes a potential gain of at most 12w, and cases 2-6 cause a potential loss of 
at least — 6(|_^J — w — logUj). These potential changes are in addition to the gain of at 
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most 331ognj — 18w discussed earlier. Thus summing the actual cost of the extract-min 
operation, c — 1 with the maximum potential gain yields an amortized cost of 89 + 391ognj. 
Summing the actual cost and the change in potential yields the amortized cost: + $j — 
< 89 + 391ogni = O {log rii). 

□ 



3.7 Global loss 

Lemma 7. - >0 

Proof. Since the structure is initially empty, $o = 0. By definition, after the execution of all 
m operations, there are no white nodes left in the collection of heaps. Since the only possible 
negative component of the potential function, the heap potential, requires white nodes to be 
negative, > and the lemma holds. □ 

3.8 Change in rank potential 

In this subsection we examine the change in rank potential in an extract-min operation. 



which was postponed in section 3.6 The proof of this bound follows in spirit the much 
simpler related proof in the original pairing heap paper [lOj, but is adjusted to take into 
account the presence of node coloring. 

Lemma 8. The rank potential gain in the EXTRACT-MIN operation where there are w 
white-white pairings is at most 181ognj — 18w 

Proof. Let s{x) be the number of white nodes in the induced subtree of x in the binary 
representation. Recall the the rank potential of a white node x is r(x) = 91ogs(a;), while the 
rank potential of a black node is 0. In this section we will also allow r(x) to be defined for 
black nodes, although it is not used in the potential function, it eases the math. 

In the original analysis of pairing heaps, using the terminology here, all nodes are white, 
and the potential function consisted only of the rank potential (without the constant of 9). In 
order to prove the O(logn) amortized runtime of EXTRACT-MIN in the original analysis, the 
change in rank potential where there are c children of the root, is shown to be 0(log?2 — c). In 
this section, this result is modified to take into account the possible presence of black nodes. 
The goal here is to show that the potential gain is 0(logn — w), where n is the number of 



white nodes and w is the number of white- white pairings in the first pairing pass. In 20 



simplified analysis of splay trees is presented, and some of our manipulations of rank potential 



are derived from this analysis, as it is easier to follow than some of the arguments in 10 . 

We first look at the analysis of a single pairing, and then proceed to analyzing separately 
the first and second paring passes. For the purpose of the analysis, we assume the first paring 
pass is done in three stages: first the white-white pairings, then the white-black pairings, and 
finally the black-black pairings. As none of these operations affect the others, this assumption 
can be made. 
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Figure 2: Two white-white pairings, one of nodes a and 6, and another of nodes d and e. 
There are no white-while pairings between h and (i, and c could be d. The node j is the first 
node in the next white-white pairing, or the rightmost white node if there is no such pairing. 

The first paring pass: White white pairings In Figure [2[ two white-white pairings are 
illustrated. Please refer to this diagram for the node naming conventions. The potential 
change of two white-white pairings is: 
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i=l 
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i=l 

LW2J-1 

= J] (2r(a,)-2r(a,+i) + 2) + 21ogn 
1=1 

<2r{ai) +w + 2\ogn 
<4 log n + ty 

The first pairing pass: Black-black pairings. These nodes have no rank potential, and 
thus nothing changes. 

The first pairing pass: Black-white pairings. There are four cases: 

Case 1: The white node starts in position a and ends in position a'. There is no 
potential change. 

Case 2: The white node starts in position a and ends in position b'. There is no 
potential gain. 

Case 3: The white node starts in position b and ends in position a'. The potential 
change for one pairing is r(a') — r{b), which summed up over k occurrences gives: 



i=l 
k-1 



<E(^(°^)-^(«m))+^(4)-^(4-i) 



i=l 



<r{a[) 
< logn 
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log r(a') — log r(n) = logr(a) — logr(6). This is at most logn for all Case 3 pairings. 
Case 4: The white node starts in position b and ends in position b'. 

This gives a potential change of logr(6) — logr(6'), which is the same potential change if 
this was a white-white pairing. Thus, this case can be treated as a white-white pairing 
for the analysis, and be incorporated into the white- white analysis presented earlier at 
no additional cost. 

Second pairing pass. Number the nodes from the top to the bottom along the right spline: 
ai, . . . Ofc. Let a'j be the node to the right of aj_i after the pairing. No matter what the 
combination of black or white nodes, or who wind the paring, the pairing involving Oj 
and a'j_^ in the second pairing pass has a rank potential change of at most log ra'j — a'-_-^^. 
Summing over all the pairing performed gives an alternating sum that is at most logn. 

□ 
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